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Abstract 

We systematically construct realistic mass matrices for the type-I seesaw mechanism out 
of more than 20 trillion possibilities. We use only very generic assumptions from extended 
quark- lepton complementarity, i.e., the leptonic mixing angles between flavor and mass 
eigenstates are either maximal, or parameterized by a single small quantity e that is of the 
order of the Cabibbo angle e ~ 9q. The small quantity e also describes all fermion mass 
hierarchies. We show that special cases often considered in the literature, such as having 
a symmetric Dirac mass matrix or small mixing among charged leptons, constitute only 
a tiny fraction of our possibilities. Moreover, we find that in most cases the spectrum of 
right-handed neutrino masses is only mildly hierarchical. As a result, we provide for the 
charged leptons and neutrinos a selected list of 1 981 qualitatively different Yukawa coupling 
matrices (or textures) that are parameterized by the Cabibbo angle and allow for a perfect 
fit to current data. In addition, we also briefly show how the textures could be generated 
in explicit models from flavor symmetries. 
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1 Introduction 



The impressive experimental advances that have been made during the past decade in 
solar [1,2], atmospheric [3], reactor [4,5], and accelerator [6] neutrino oscillation experiments, 
have very well established that neutrinos are massive. Since neutrinos are massless in the 
Standard Model (SM), the observation of neutrino masses provides evidence for new physics, 
such as an underlying Grand Unified Theory (GUT) [7] (see also Ref. [8]). It is therefore 
believed that the smallness of the absolute neutrino mass scale m v ~ 10~ 2 . . . 10 _1 eV 
compared to the electroweak scale ~ 10 2 GeV gives us important information on the nature 
of the new physics. Today, the most widely accepted mechanism to generate small neutrino 
masses is the seesaw mechanism [9, 10], in which the smallness of neutrino masses is linked 
to the hierarchy between the electroweak and the GUT scale Mqut — 2 x 10 16 GeV [11]. 

In the type-I seesaw mechanism [9], the set of SM neutrinos v\ {% — 1, 2, 3 is the generation 
index) is extended by three right-handed neutrinos uf, which are total singlets under the 
SM gauge group Gsm = SU(3) C x SU(2) L x U(l)y- In the basis (z/i, v 2 , v\i v 2i v t)i this 
leads after electroweak symmetry breaking to a complex symmetric 6x6 matrix 



where 0, Mr>, and Mr are 3x3 matrices. The upper left matrix has zero entries since there 
is no Higgs triplet that could have directly coupled to the V{. The entries in Mb are protected 
by electroweak gauge invariance and they are therefore of the order ~ 10 2 GeV, while the 
matrix elements of Mr are of the order of the B — L breaking scale Mb-l — 10 14 GeV. 
After integrating out the right-handed neutrinos, we arrive at the effective low-energy 3x3 
neutrino Majorana mass matrix 



which gives rise to neutrino masses of the order m v ~ 10 2 eV. The seesaw mechanism is 
attractive because Mb-l is very close to Mqut, indicating a GUT-origin of neutrino masses. 

In GUT models, quarks and leptons are unified into multiplets, which is known as quark- 
lepton unification and one possibility to explore GUTs at present energies is to search for 
signatures of quark-lepton unification in the fermion mass and mixing parameters. Most 
notably, quark-lepton unification has to give an answer to the question why the quark mixing 
angles in the Cabibbo-Kobayashi-Maskawa (CKM) matrix Vckm [12] and the leptonic mixing 
angles in the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix £7pmns [13] are strikingly 
different. In the quark sector, all CKM mixing angles are small and can be approximately 
written as powers of the Cabibbo angle Qq. In contrast to this, in the lepton sector, only 
the reactor angle #13 is small, whereas the solar angle Qyi and the atmospheric angle #23 
are both large. In addition, while the quark and charged lepton mass ratios are strongly 
hierarchical, the neutrino masses exhibit, if any, only a mild hierarchy (for a recent global 
fit of neutrino data see, e.g., Ref. [14]). 

Recently, quark-lepton complementarity (QLC) [15] (for an early approach see Ref. [16]) has 
been proposed as a possibility to account for the differences between the quark and lepton 
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mixings. In QLC, the quark and lepton mixing angles are connected by the QLC relations 

012 + C « tt/4, 23 + cfe « tt/4, (3) 

where 8 c t> = arcsin V c b- The crucial observation is that sum rules of the types shown in Eq. (j3J) 
can be easily obtained when the mixing among the neutrinos and among the charged lep- 
tons is described by maximal or CKM-like mixing angles. In this way, for example, the 
observed value of the solar angle 9 12 ~ 33° could be understood in terms of maximal (7r/4) 
and Cabibbo-like (0c) mixing in the individual neutrino and charged lepton sectors. A 
complementary approach to the solar angle seems, on the other hand, to be suggested by 
the tri-bimaximal mixing scheme [17]. The properties of QLC have been studied in various 
respects: as a result of deviations from bimaximal mixing [18], in connection with sum 
rules [19], with emphasis on phenomenological implications [20], together with parameter- 
izations of C/pmns m terms of 0c [21], in view of statistical arguments [22], in conjunction 
with renormalization group effects [23], and in model building realizations [24]. 

In Ref. [25], we have proposed an extended QLC, in which the mixing angles in both the 
charged lepton and the neutrino sector can take any value in the sequence 7r/4, e, e 2 , . . . , 
where e is of the order the Cabibbo angle e ~ 9c- In this paper, we will implement extended 
QLC in the type-I seesaw mechanism by assuming that all the mixing angles of charged 
leptons and left- and right-handed neutrinos take their values in this sequence. We also 
suggest that the mass eigenvalues of and Mr are described by powers of e as well. In this 
approach, the observed large mixing angles 9 12 and 023 can come from the charged leptons 
and/or neutrinos]]] Moreover, in the neutrino sector, large mixing angles can originate from 
Mr, and/or Mr. We systematically search for all mass matrices of charged leptons and 
neutrinos that satisfy the extended QLC assumptions, and extract all solutions that are 
consistent with current data in the CP conserving case. 

The paper is organized as follows: We first motivate the assumptions underlying the ex- 
tended QLC approach from the phenomenological and model building point of view in 
Sec. This section can be skipped by the reader already familiar with extended QLC. 
Next, in Sec. [31 we describe the method for constructing all valid charged lepton and seesaw 
mass matrices that are compatible with data, demonstrate how to obtain the corresponding 
textures, and address further properties of our procedure. For the normal neutrino mass 
hierarchy, we first discuss the full sample of all valid mass matrices in Sec. HI and then we 
show a selection of textures in Sec. A qualitative discussion of the inverted and degen- 
erate neutrino mass schemes is included in Sec. [6] and a summary and conclusions can be 
found in Sec. [7J Details of our method can be found in Appendix [A] and [B] 

2 Motivation 

In this section, we first present a brief review of the observed hierarchies of fermion masses 
and mixing angles and relate them to a small expansion parameter e that is of the order of 
the Cabibbo angle. Then, we discuss two representative SU (5) GUT examples that obtain 

1 For a recent study on the reconstruction of the seesaw mechanism from low-energy data see Ref. [26] 
and large mixing angles coming from the charged lepton sector were also considered, e.g., in Ref. [27]. 
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the observed mass and mixing parameters from flavor symmetries. The observations made 
here will later, in Sec. 13.21 serve as a motivation for the hypotheses in extended QLC. 

2.1 Masses and Mixings of Quarks and Leptons 

One of the most striking features of the fermion sector is that the mass and mixing param- 
eters of quarks and charged leptons are strongly hierarchical. It is well known that these 
hierarchies can be approximately described by a small number e ~ 0.2. In the Wolfenstein 
parameterization [28], for example, the CKM matrix is given by 

/ l-\e 2 e A(p-i V )e 3 \ 

Vckm = ~e l-\e 2 Ae 2 , (4) 

\A(1 -p-ir])e 3 -Ae 2 1 / 

where e is of the order of the Cabibbo angle 9q ~ 0.2, and A, p, and t], are order unity 
parameters (for an update see Ref. [29]). Order of magnitude wise, the quark mixing angles 
can be written in terms of the parameter e as 

|K,|~e, \V cb \~e 2 , \V ub \~e 3 . (5) 

An interesting connection between the Cabibbo angle and the quark masses is established 
by the Gatto-Sartori-Tonin-Oakes relation 6q = \fm^Jrn s [30], suggesting that also the 
fermion mass ratios arise from powers of e. In fact, the mass ratios of the up and down 
quarks can, e.g., be crudely represented as powers of e as^| 

m„ : m c : m t = e 6 : e 4 : 1, : m s : m b = e 4 : e 2 : 1, (6) 

where m b : m t ~ e 2 , m T : ~ 1, and m t ~ 175 GeV, whereas the mass ratios of the 
charged leptons are crudely given by 

m e : m M : m r = e 4 : e 2 : 1. (7) 

These mass ratios have all to be understood as order of magnitude relations and depend 
on the energy scale. At one loop, in the minimal supersymmetric standard model (MSSM), 
the only changes in these relations are m c : m t ~ e 4 and rrib '■ rn t ~ e 3 at ~ 10 15 GeV (see, 
e.g., Ref. [31]). The changes due to renormalization group (RG) running are thus only 
comparatively small. 

In the neutrino sector, we have a situation that is substantially different from the charged 
fermion sectors. Experimentally, the PMNS matrix reads (c/. also Ref. [32]) 

/0.82 - 0.85 0.52 - 0.57 < 0.12 \ 
t/pMNS = 0.26 - 0.49 0.48 - 0.68 0.65 - 0.76 , (8) 
\0.27- 0.49 0.49 - 0.69 0.64- 0.75/ 

which has, unlike Vckmj large off-diagonal entries. In the standard parameterization, the 
la ranges for the solar and the atmospheric mixing angles are then [14] 

sin 2 fl 12 = 0.30^, sin 2 # 23 = 0.5^°;°?, (9) 
2 We are interested here in an SU(5) compatible fit. 
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whereas we have a 3<j upper bound on the reactor angle that is sin 2 #i3 < 0.041. The best 
fit values of the mixing angles correspond to maximal atmospheric mixing # 2 3 ~ vr/4 and 
large, but not maximal, solar mixing # 12 ~ vr/4 — 0q, and a small reactor angle 13 < 0q. 
The la bounds on the solar and atmospheric mass squared differences are [14] 

Ami = (7-9^:1) x KT 5 eV 2 , Am 2 tm = (2.5^; 2 5 ) x lO^eV 2 . (10) 

The sign of Ami is positive, whereas the sign of Am 2 tm can be either positive or negative, 
leading to currently three possible types of allowed neutrino mass spectra. Note that we 
roughly have 

Ami : Aml tm ~ e 2 . (H) 

It is thus plausible that the neutrino sector, just like the quarks and charged leptons, 
is described by the same control parameter e. The order of magnitude relations for the 
neutrino masses may thus be written as 

mi : m 2 : m 3 = e 2 : e : 1, mi : m 2 : m 3 = 1 : 1 : e, mi : m 2 : m 3 = 1:1:1, (12) 

where m 1; m 2 , and m 3 , denote the masses of the 1st, 2nd, and 3rd neutrino mass eigenstate. 
The 1st, 2nd, and 3rd equation in Eq. ffl2l) describe a normal hierarchical (NH), inverse 
hierarchical (IH), and quasi degenerate (QD) spectrum, respectively^ 

We thus see that the CKM angles and mass ratios of quarks and leptons are roughly given 
by some power ~ e n of the Cabibbo angle e ~ 9c- Besides that, the phenomeno logical 
QLC relations 623 ~ tt/4 — e 2 and 812 ~ 7r/4 — e involve maximal mixing angles. In the 
next section, we will dicuss how these mass and mixing parameters may be reproduced in 
models. 

2.2 Examples with Quark-Lepton Unification 

The common appearance of the masses and mixing angles in the quark and lepton sectors 
may point towards a quark-lepton unified theory. One might therefore wonder whether a 
description of the fermion mass and mixing parameters as given in Sec. 12.11 can indeed be 
obtained in explicit models. We are interested in models in which the fermion mixing angles 
- prior to going to the mass eigenbasis - can be maximal or are given by some power of the 
Cabibbo angle e ~ 9c which also describes all fermion mass ratios. 

For this purpose, let us briefly review two GUT models [33,34] based on 577(5), which 
show that in a quark-lepton unified theory it is (i) actually possible to generate realistic 
hierarchical fermion mass ratios and mixing angles that are described by powers ~ e n and 
that (ii) one can predict exact maximal mixing compatible with these hierarchies. We view 
the two GUT models as two representatives of a broad class of possible realistic models using 
Abelian (for early work see Ref. [35] and for more recent models see, e.g., Refs. [31,36]) or 
discrete non- Abelian (for recent studies including the quark sector see, e.g., Ref. [37] and 
in the context of GUTs see, e.g., Ref. [38]) flavor symmetries^ Our observations can be 
viewed as a further motivation for the definition of extended QLC in Sec. 13.21 where we 

3 For NH neutrinos, one can compute e from the current best-fit values, which gives 0.15 < e < 0.22 (3a). 
4 For a more complete list of references on discrete non- Abelian flavor symmetries see Ref. [39]. 
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Figure 1: Generation of higher-dimension fcrmion mass terms via the Froggatt-Nielsen mechanism. 
ipL and ipn are left- and right-handed SM fcrmions. Internal solid lines denote superheavy fermions with 
common mass Mp, the circled cross is the usual SM Higgs vacuum expectation value (VEV), whereas the 
crosses without circle represent universal VEVs v of SM singlet scalars that break the flavor symmetry. 
After integrating out the heavy fermions, the resulting mass term is then given by the effective dimension-n 
operator (H)e n ip L il)R, where e = v/Mp serves as a small expansion parameter. 



will, in particular, claim that both the neutrinos as well as the charged leptons can exhibit 
mixing angles that are maximal and/or ~ e n . In this way, the size of 612 will then be simply 
the result of a QLC-type sum rule. 

Example 1: Cabibbo-type mass and mixing hierarchies Our first example is a 
super symmetric SU(5)' x SU(5)" GUT with a U(1) N = 11^11(1) j flavor symmetry group 
[33]. It yields, as a result of U(1) N flavor symmetry breaking, the masses and mixing angles 
of quarks and leptons roughly as powers of e. The ith generation (i = 1, 2, 3) is charged under 
SU (5)' x U (l)ki labels a suitable subgroup U(l)k t C U^ =1 U(l)j) in an SO (10) compatible 
way as 10(— l)j + 5(3)j+l(— 5)j, where the numbers in parenthesis denote the 11(1)^ charges 
of the respective multiplets and the 577(5) singlets are the right-handed neutrinos. Higgs 
superfields break the U(1) N gauge symmetry such that masses for quarks and leptons arise 
from higher-dimension terms via the Froggatt-Nielsen mechanism [40] illustrated in Fig. [TJ 
In the presence of SM singlet scalar "flavons" , that break the flavor symmetry by acquiring 
universal vacuum expectation values (VEVs) v (crosses), and superheavy fermions with 
common mass Mp, which are charged under the flavor symmetry (internal solid lines), 
the mass terms of quarks and leptons become suppressed by integer powers of a small 
parameter e = v/Mp that controls the flavor symmetry breaking. The integer power of e 
is solely determined by the quantum numbers of the left- and right-handed fermions ipL 
and ipR under the flavor symmetry. As a result, the mass matrices of the up quarks, down 
quarks, charged leptons, and neutrinos become 

e 5 e 3 \ /e 4 e 6 e 3 \ /e 4 e 2 e 

M u = m t I e> e 4 e 2 , M d = Mj = m b e w e 2 e 2 , M eS = m A e 2 1 1 
e 2 1/ \e 8 e 3 1/ \ e 1 1, 

(13) 

The representation of the mass matrices in Eq. ([TBI are examples of what we will call in the 
following textures: These are descriptions of the mass matrices showing only the order of 
magnitude (up to order one Yukawa couplings) of the entries in terms of powers of a small 
number e that parameterizes the flavor symmetry breaking. At an order of magnitude level, 
the textures in Eqs. (fT3l predict for the quarks and charged leptons the mass ratios and 
mixing angles of Eqs. (jSJ), (EJ), and whereas the neutrino mass spectrum is of the type 
mi : m>2 : = e : e : 1. The reactor mixing angle is small and of the order #13 ~ e, while 
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the solar and the atmospheric mixing angles are large and of the orders #12 ~ 1 and #23 ~ 1. 
A proper choice of the order one Yukawa couplings then allows to reproduce the solar and 
atmospheric mixing angles close to the current best fit values. 

Example 2: Maximal mixing - Our second example is an SU(5) GUT with a non- 
Abelian discrete flavor symmetry between the 2nd and 3rd generation [34]. It predicts a 
maximal atmospheric mixing angle #23 = vr/4 and a zero reactor angle # 13 = as a conse- 
quence of the discrete symmetry. The basic flavor symmetry of the model is a discrete Z 2 
exchange symmetry that implements a maximal atmospheric mixing angle and acts on the 
SU(5) multiplets of the 2nd and 3rd generation as Z 2 : 52 <-» 53, IO2 <-» IO3, I2 <-> I3, 
where the subscripts denote the generation indices. In addition, the model has a U(l) family 
number symmetry that does not commute with the above Z2 generator. As a consequence, 
the resulting down quark and charged lepton mass matrices can accommodate the hierarchi- 
cal down quark and charged lepton masses and the CKM and PMNS mixing angles arise in 
the up quark and the neutrino sector, respectively. The total non-Abelian flavor symmetry 
enforces in M c g a fi — r exchange symmetry. This predicts a maximal atmospheric mixing 
angle #23 = j and a vanishing reactor angle #13 = 0. The solar angle 612, on the other hand, 
is large and can easily reproduce the current best fit value. 

These two examples show that hierarchical masses and mixings described by powers of e as 
well as maximal mixing can be predicted in explicit GUT models from flavor symmetries. As 
a result of this motivation section, it is therefore plausible to assume that all mixing angles 
and mass hierarchies in the quark and lepton sectors are generated by a single small quantity 
e ~ 8q augmented by possibly maximal mixing in the up- and/or down-type sectors. This 
basic observation will be the basis for our hypotheses underlying extended QLC. 

3 Method 

In this section, we introduce our method for implementing extended QLC in the seesaw 
mechanism. For this purpose, we first briefly review the seesaw mechanism and discuss our 
notation for parameterizing the mass and mixing parameters in Sec. 13.11 Next, we define 
our QLC assumptions in Sec. 13.21 and outline our approach for generating and selecting the 
mass matrices of charged leptons and neutrinos in Sec. 13.31 i.e., we describe our general 
procedure. While Sec. 13.31 is somewhat qualitative in some points, we give more details and 
a comment on the complexity in Appendix [A] As the next step, we demonstrate how we 
produce textures in Sec. 13.41 Finally, we discuss the role of mass ratios, input parameters, 
and RG running in our routine in Sec. 13.51 

3.1 Neutrino Mass and Mixing Nomenclature 

In what follows, we assume that the left-handed SM neutrinos acquire their masses via 
the type-I seesaw mechanism [9]Jj In the type-I seesaw mechanism, the Yukawa couplings 

5 The type-II seesaw mechanism [10] in extended QLC has already been covered in a previous analysis 
of the possible effective 3x3 Majorana neutrino mass matrices in Ref. [25] that can be viewed as being 
generated by the coupling to some Higgs triplet with small (~ 1CT 2 cV) VEV. 
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generating the charged lepton and neutrino masses are 

C Y = -(Y&jH*^ - (Yo^m^ - X -[M R ) i3 vlv] + h.c., (14a) 

where = (z^, ei) T , e\, and uf (i = 1,2,3 is the generation index), are the left-handed 
lepton doublets (£«), the right-handed charged leptons (ef), and the right-handed SM singlet 
neutrinos (vf). In Eq. (I14al) . H is the Higgs doublet, Yi and Yd are the 3x3 Dirac Yukawa 
coupling matrices of the charged leptons (Yg) and neutrinos (Y^), M R is the 3x3 Majorana 
mass matrix of the right-handed neutrinos, and icx 2 is the 2x2 antisymmetric tensor. After 
electroweak symmetry breaking, H develops a vacuum expectation value (H) = (0, v/\/2) T , 
where v ~ 10 2 GeV, and the mass terms of the leptons become 

£ma SS = -{M^e] - (Md)^ - -{M R ) i3 v?v? + h.c, (14b) 

where Mi = {H)Yi is the charged lepton and Mp = (H)Y^ the Dirac neutrino mass matrix. 
Mi and Mp are complex 3x3 matrices that are described by 18 parameters and have 
entries of the order ~ 10 2 GeV. The matrix Mr is complex, symmetric, and described by 12 
parameters. It has matrix elements of the order the B — L breaking scale Mb-l ~ 10 15 GeV. 
The resulting complex symmetric 6x6 neutrino mass matrix is given in Eq. (Tj[|), and after 
integrating out the right-handed neutrinos, this gives the effective 3x3 neutrino mass matrix 
M c fj in Eq. ([2]) leading to masses ~ 10~ 2 eV for the active neutrinos. 

To analyze the origin of leptonic mixing, we consider the diagonalization of the mass terms 
in £ mass in Eq. (114bp by unitary matrices. Using the convention in Ref. [25], we can always 
write a general unitary 3x3 matrix f/ un i tary as 

Unitary = diag (e^ 1 , , ) ■ U ■ diag (e ia \ e ia \ 1) , (15a) 

where the phases <pi, ip 2 , <f3, «i, and a 2 , take their values in the interval [0, 27r] and 

' c 12 c 13 s 12 c 13 si 3 e" 1? \ 

U = 



-S12C23 - c 12 s 2 3Si 3 e^ c 12 c 23 - s 12 s 23 s 13 e lS ^ s 23 c 13 
\ s 12 s 23 - c u c 23 s 13 e lS -c u s 23 - s 12 c 23 s 13 e lS c 23 c 13 ) 



(15b) 



is a CKM-like matrix in the standard parameterization with sy = sin 9ij, = cos 9^, where 
Oij e {9 12, 9\ 3 , 9 23 } lie all in the first quadrant, i.e., 9^ e [0, |] , and 5 £ [0, 2tt]. The matrix 
U is thus described by 3 mixing angles 9j~ and one phase 5, i.e., it has 4 parameters. The 
matrix t/ un i t ary has five additional phasesj and contains therefore in total 9 parameters. 

The leptonic Dirac mass matrices Mi and Mjj, and the Majorana mass matrices Mr and 
M e g are diagonalized by 

Mi = UiMf^Ul M D = U D M d ^U j D „ M R = U R M^U T R , MeB = U V M%*U? , (16) 

where Ui, Ui>, Up, Ue>>, U r , and U u , are unitary mixing matrices, whereas M^ iag , M D iag , M^ iag , 
and M^ g , are diagonal mass matrices with positive entries. We can always write the mixing 
matrices as the products 

U x = D X U X K X , (17) 



'For recent discussions of rephasing invariants in the lepton sector see, e.g., Ref. [41]. 
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where U x are CKM-like matrices that are parameterized as in Eq. (115bj) , while D x and K x are 



given by D x = diag(e lip i , e 1 ^, e 1 ^) and K x = diag(e 1Q i , e ia % , 1), where the index x runs over 
x = £,£',D, D', R, v. The phases in D x and K x are all in the range <pf, tp^ , ip^ , af , a\ G [0, 27r] . 
Each of the matrices U x in Eq. (fT7j) contains four mixing parameters: three mixing angles 
and one phase. We denote the parameters of U x by ^12 ? ^13 ? ^23 5 an d & x - For each of the 
matrices U x in Eq. (fT7j) . we define the mixing parameters by identifying in Eq. (115b|) the 
mixing angles as 6^ — > 9fj, and the phase as 5 — > £ a . 

The PMNS matrix is given by 

t^PMNS = C^C/v = f^PMNS-^Maj, (18) 



where U is a CKM-like matrix parameterized as in Eq. (115b[) . and i^Maj = diag(e 1 ^ 1 , e 1 ^ 2 , 1 



contains the Majorana phases <p\ and (f>2- The CKM-like matrix t/pMNS m Eq. (118j) is 
described by the solar angle 612, the reactor angle #13, the atmospheric angle #23, and the 
Dirac CP-phase 5, which we identify in the standard parameterization of Eq. (115bj) as 



$ij — > 9ij and 5^5. The PMNS matrix has thus 3 mixing angles and 3 phases and contains 
therefore 6 physical parameters. 

Let us next express M e g- in terms of the mass eigenvalues and mixing angles introduced 
above. Inserting Eq. (ITT)) into Eq. (1TE1) . we find 

M D = D D U D K D M% as K* D ,ff D ,D* D ,, (19a) 

M~ R l = D^KUM^T'KrU^. (19b) 

The effective neutrino mass matrix M e g in Eq. (j2J) can thus be written as 

= —D D U D KM^ g Jj] D ,DU^ i (i^) 2 (Af^ iag ) _1 f/^^C/J, Mp^KU^D D , (20a) 

where we have introduced K = K* d Kdi and D = D* D ,D* R . We have denoted the parameter- 
ization of M e fj in Eq. (I20al) by an extra superscript "th" for "theoretical" , since none of the 
mass and mixing parameters on the right-hand side of Eq. (|20aj) are directly measurable in 
neutrino oscillations. Note that in the CP conserving case, the matrix {K^) 2 drops out of 
the expression for M e s in Eq. (120aj) . Equivalently to Eq. (I20aj) . using Eqs. (ITS")) and (ITTj) in 
the expression for M e g in Eq. (TIE)) , we can write M e g- also in the parameterization 

M e 7 = D e U £ K £ U PU mK^M^ s U^ um K e UjD e , (20b) 

where we have chosen the superscript "exp" for "experimental" , to label the representation 
of M cS in Eq. (I20bl) . since M^ p involves the matrices M^ g and Upuns containing the 
experimentally accessible mass and mixing parameters. Note that in the CP conserving 
case, (i^Maj) 2 drops out of the expression for M^ p . It is clear that = M^ v = M e g, 
since and are just different parameterizations of M e fj. 

3.2 Hypotheses for Extended QLC 

We will now formulate the assumptions underlying extended QLC in the type-I seesaw 
mechanism. Motivated by the discussion of the GUT examples in Sec. 12.21 extended QLC 
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will in this context include assumptions on (i) the mixing parameters of Ui,Ur,Uri, and 
Ur, and assumptions on (ii) the eigenvalues of the mass matrices Mi, M e g, Mr, and Mr. 

Mixing angles - Consider first the possible mixing angles in extended QLC Following the 
idea of Ref. [25], we begin by assuming that all mixing angles parameterizing in Eq. (fTol) 
the mixing matrices Ue, Ur>,Ur,i, and Ur, can a priori take any of the values in the sequence 
7r/4, e, e 2 , . . . , where e is a small number. Motivated by the quark sector, we will take 
e ~ 0.2, i.e., we assume that e is of the order the Cabibbo angle. This applies the concept 
of extended quark- lepton complementarity in Ref. [25] to the mixing matrices Ud, Uri, and 
Ur, that diagonalize the renormalizable neutrino mass terms in Eq. (114b|) . Since the current 
la error on the leptonic mixing angles is at most of the order e 2 (see, e.g., Ref. [14]), we 
will truncate the sequence of mixing angles 7r/4, e, e 2 , . . . after the element e 2 and identify 
there all terms of the order e n with n > 3 simply by "0" . In other words, in extended QLC, 
we restrict the possible range of mixing angles 9fj (x = £, D, D' , R) to the set of values 
9fj G {7r/4, e, e 2 , 0}, where "0" represents mixing angles ~ e n with n > 3. Since we want to 
compare with current neutrino data, the assumptions on the mixing angles are formulated 
at low energies ~ 1 GeV. The impact of RG running when assuming these angles at a high 
scale will be discussed later in Sec. 13.51 Note that we cannot just rotate away Ue, because 
this would make our mixing angle assumptions as powers of e meaningless if induced by an 
underlying theory, such flavor symmetry. 

Mass eigenvalues - Next, let us specify the types of mass eigenvalues that we assume in 
extended QLC. We have to distinguish two types of mass eigenvalues - those of the lepton 
mass matrices and M e g in the low-energy effective theory, and those of the Dirac and 
heavy Majorana mass matrices Mr, and Mr of the neutrinos. 

We assume that the mass spectra in the low energy effective theory are those of Eq. ([7j) and 
Eq. (I12p . That means we assume m e : m^ : m T = e 4 : e 2 : 1 for the charged leptons, whereas 
we have in the neutrino sector mi : rri2 : = e 2 : e : 1 for a NH, mi : rri2 : = 1 : 1 : e for 
an IH, and mi : m 2 : m 3 = 1 : 1 : 1 for a QD neutrino mass spectrum. As we will discuss 
later, our results are actually completely independent from the details of the charged lepton 
spectrum. 

Let us now consider the mass eigenvalues of the Dirac and Majorana mass matrices Mr, and 
Mr. We denote the mass eigenvalues of Mr, by mf, , and mf, and the mass eigenvalues 
of Mr by mf,m^, and mf. Notice that these are not directly observable at low energies, 
but as before, we write the mass eigenvalues of Mr, and Mr as powers of e. Similar to 
Eqs. (|7|) and (jl2|) . we will therefore parameterize the mass eigenvalues Mr, and Mr as 

mf : mf : mf = e a : e b : e c and mf : mf : mf = e a ' : e b ' : e c ', (21) 

where a, b, c, a 1 , b', and c', are suitable non-negative integers < 2, and we define the absolute 
mass scales by = rriR,e c and m^ = Mr_r€ c . As we will discuss in Appendix IA.1[ one 
can restrict the possible range of the integers in Eq. (j2ip to a fairly small set of numbers, 
such that it is sufficient to test all possible combinations up to second order in e. 

In total, we see that in our hypotheses for extended QLC all mass hierarchies and small 
mixing angles become described by powers of the hierarchy parameter e. In Sec. 13.51 we will 
discuss the validity and precision of our assumptions in extended QLC when taking, e.g., 
GUT relations among fermion masses and RG effects into account. 



9 



3.3 Procedure: Systematic Construction of the Parameter Space 



„ , , . v oscillation and mass t. «- , , Flavor 

Observation: , Model: symmetry, 

experiments, ... GUT. 
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Seesaw realizations: M;, Mjj, Mr 



Reduction 



Textures 



O 
Cf 

£ 

0) 




Figure 2: Procedure for obtaining the seesaw realizations and texture sets in extended QLC. 



Let us now describe our three-step procedure for generating all textures Mg, Mp, and M^, 
which satisfy extended QLC for the seesaw mechanism (cf., Fig. [2] for illustration): 



First step - We generate all possibilities for the effective neutrino mass matrix M*^ in 
Eq. (I20ap . Here, we assume that the mixing angles entering Up, Uc, and Ur, can take any 
values in the set 

0*. eW4,e,e 2 ,0}, (22) 

where x = D, D', R. Moreover, we suppose in Eq. ([20aD that M^ iag and M^ iag are on the gen- 
eral forms as given in Eq. ( 12"TT) with eigenvalues 1, e ore 2 . For simplicity, we will confine our- 
selves to the CP conserving casd^ were all phases are taken from the set 5 X , ipf, y?f , ^3, a\, af G 

{0,7T}. 

Second step - We generate all possibilities for the neutrino mass matrix M^ p in Eq. fl20b[) . 
For £/pMNS; we use values motivated by the current best-fit values. In particular, we use the 



7 Note that in the case of CP violation, some textures may change due to cancellations, so the number 
of textures will increase, like one would expect. Nevertheless, a complete systematic analysis (all phases 
between and 2tt are allowed) is up to now not possible because of lack of computing power. This may 
change in some years. 
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following input, which could be experimentally confirmed or rejected in the coming yean 




choice for these parameters can be equally well applied, but it will change the final results. 

For Ug, we follow Eq. (1221) . and for the phases, we test all real possibilities. Furthermore, we 
insert the neutrino mass spectra given in Eq. ([T2D into M^ e , and again test all possibilities. 

Third step - Next, we match all possibilities from step 1 and step 2, i.e., we select all 
parameter combinations for which 



at 0(e 3 ). For details (and the exact numerical implementation) of the matching procedure, 
see Appendix IA.21 in particular, note that our procedure automatically factors out the 
overall neutrino mass scale, i.e., m v = m 2 D /Mn will be automatically satisfied. In the 
following, we will call a (seesaw) realization a valid set {Md,Mr, Ug}, or more precisely, 
a combination of all involved mixing parameters, phases, and mass hierarchies, for which 
Eq. (1241) is fulfilled. In other words, a realization is thus a set of input parameters compatible 
with current experimental best-fit values which describes and the 6x6 matrix M v 
completely, and it contains the left-handed charged lepton mixing. In total, our procedure 
requires that we systematically scan 20 trillion different possible realizations. For details on 
the complexity, see Appendix IA.3I 

Note that in our procedure the observed large leptonic mixing angles can be generated either 
in the charged lepton sector and/or the neutrino sector. Furthermore, in the neutrino sector, 
large mixings can arise from the Dirac neutrino and/or the Majorana mass matrix of the 
right-handed neutrinos. This means that we not make any special assumptions simplifying 
the structure of the seesaw mechanism, such as taking Mg or Mr to be diagonal, or Mr, 
symmetric. 

3.4 Texture Extraction and Order Unity Couplings 

Let us now describe how we extract the textures for the charged leptons and neutrinos 
from the seesaw realizations that satisfy Eq. (|2~4"1) . In the course of applying the three step 
procedure described in Sec. 13.3} we have already produced for each realization the pair of 
matrices Mo and Mr. In the same way, we determine for each valid realization the charged 
lepton mass matrix 10 ! Mi by rotating to the left-handed flavor basis as Mg = UgM^ iag , where 
M^ iag contains the masses given in Eq. ([7]). Next, we analytically expand the mass matrices 
in e as 



8 See, e.g., Refs. [42,43] for long-baseline experiments on a scale of the coming ten years, Refs. [44,45] for 
an up scale reactor #12 measurement, Rcf. [46] for the potential of various different supcrbeam upgrades, 
and Rcf. [47] for a neutrino factory measurement. 

9 Note that these values correspond not only to the best-fit values, but are also often considered as an 
interesting symmetry limit in "exceptional" [39] neutrino mass models. 

10 Notc that we choose, for simplicity, the right-handed charged lepton mixing matrix to be the unit matrix 
Ui> = 1. This choice is, however, not limiting our procedure since Ug> docs not enter into E/pmns- 
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(24) 
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Figure 3: Distribution of the order one couplings in Mi (left), Mp (center) and Mr (right) for all valid 
seesaw realizations of NH neutrinos (of all orders and for all matrix elements). The gray-shaded region 
marks the area in which the order one couplings are smaller than e. 

where x = £,D,R, and identify for each matrix element the leading contribution as the 
lowest order in e. The texture for M x is then found by substituting each matrix element of 
M x by its leading power e k . For k > 3, we take e k — > 0. In doing so, we always "drop" the 
order one coefficients that multiply the leading powers e k and round them to one (unless they 
are zero). For a given realization, we call the collection of textures for Mi,Md, and Mr a 
texture set. For any texture set obtained in this way, we argue that one can always adjust 
the order unity coefficients such that the Yukawa couplings are brought in perfect agreement 
with data. An important prerequisite for this statement is that the different orders in the 
expansion in Eq. (I25j) do not "interfere" with each other, i.e., the involved coefficients are 
really of order unity. This property for the coefficients will be checked below. It is important 
to keep in mind that, in general, more than one realization may lead to the same texture 
set. Because of the reduction from in general several realizations to one texture set, we will 
call our texture set producing technique "texture reduction" . Note also the important fact 
that the texture reduction is based on an analytic expansion in e as opposed to just using 
powers of 0.2 for a purely numerical fit of the mass matrix elements. 

As an example for texture reduction, consider the following mass matrix Mo (cf., tex- 
ture/realization #1 in Tabled]): 



— e 



\ 



V2 2V2 V2 2V2 

2 ,2 J2 



c 2 
c 2 
c 2 




Here, "— ►" symbolizes, up to an overall mass scale, the identification of the leading order 
terms in the expansion in e that contribute to the mass matrix elements in Md- The matrix 
on the right hand side of — > then represents the texture corresponding to the mass matrix 
Md on the left hand side. 

In order for the expansion in Eq. (j2"51) to be useful, the Yukawa couplings should be of order 
unity (at least if one wants not to rely on some sort of fine-tuning). The expansion suggests 
a criterion for what "order unity" actually means: the coefficients should lie in the interval 
[e, . . . , e _1 ] ~ [0.2, . . . , 5], such that they would not imitate different orders in e. In Fig. [HJ 
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we show to which extent this condition is indeed satisfied for the valid realizations. Fig. [3] 
depicts the coefficients of the in Eq. (1251) for Mi (left), M D (center) and M R (right) 
for all valid seesaw realizations and NH neutrinos (including all orders and for all matrix 
elements). Interestingly, it turns out that in 99.9% of all cases, the order unity Yukawa 
couplings lie in the range e, . . . , e _1 , which justifies the expansion of the textures in e. For 
example, for Mi and Md, the peaks at around 0.7 ~ l/y/2 are predominant. Note that, 
though in Fig. [3] the first bin corresponds to coefficients smaller than e, our mapping of 
textures is unambiguous. This is because in the cases where the leading order coefficient of 
a matrix element becomes very small (< 0.2), the other coefficients become also very small 
- with the exception of 0.10% of all cases for the matrix element si"! As a consequence, the 
leading order term indeed rarely numerically interferes with the higher orders. This picture 
would change if e was larger. For instance, for e = 0.5, the leading order identification for 
Md fails in 0.43% of all cases, and for Mr it fails in 0.01% of all cases, which is, however, 
still only a remarkably small fraction of all cases. 

3.5 Dependence on Input and Renormalization Group Effects 

Let us now discuss how the choice of the input parameters from measurement affects our 
results. First of all, we use the current best-fit values as an input (cf., Eq. (1231) ). Any other 
choice of experimental input parameters will give different results, but our procedure can 
be equally well applied to any other preferred choice of input values. We have chosen these 
input values because these results can be expected to remain valid for the next ten years 
or so, unless the best-fit values change (such as if sin 2 2^13 ~ 0.1 were indeed found). We 
have also computed the dataset for different values of #13, but a presentation of these results 
would clearly exceed the scope of this paper. 

As far as the charged lepton mass spectra in Eq. (|7]) are concerned, a variation of the mass 
ratios will not have any effect at all on the selection and extraction of the texture sets: 
Any change in the spectrum drops completely out of our routine. The particular choice of 
the charged lepton mass hierarchies in Eq. (J7|) has been motivated by comparison with the 
down quark spectrum in SU(5). Choosing a different parameterization of the charged lepton 
masses, for example with different powers of e, will of course have an effect on the form of 
the extracted texture sets for Mi, but it does not affect the selection of the realizations in 
any way. In particular, a modification of the charged lepton mass spectra in Eq. (j7j), e.g., 
to implement the Georgi-Jarlskog relation m M : m T = 3m s : [48], would not change any 
of our results for the textures. In fact, one can easily obtain the textures of the charged 
leptons for any other choice of the hierarchy by using Ut directly. 

The extended QLC hypotheses in Sec. 13.21 should hold at high energies such as Mgtjt, but 
they are compared in our method with low-energy data. We therefore have to address the 
stability of the extended QLC assumptions under RG running from the high scale down to, 
say, around ~ 1 GeV. Generally, the empirical QLC sum rule 12 = 7r/4 — 9 C is satisfied 
up to a precision of about ~ 1°. It is therefore reasonable to take in our routine nonzero 
mixing angles 6fj into account that can be as small as e 2 ~ 2°. Moreover, it is known that 
the Cabibbo angle 8c does practically not run [49] , and V c b ~ e 2 changes typically only by a 

11 All these exceptions appear only in Md- 
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factor smaller than 2 when running from ~ 1 GeV up to the Planck scale ~ 10 19 GeV [50]. 

Let us have a more precise look at the RG evolution of neutrino masses and mixings [51]. 
First, note that, due to the smallness of the charged lepton Yukawa couplings, the running 
of a possibly maximal atmospheric mixing angle #23 is negligible, unless one works in the 
MSSM with large tan/? [43]. Simple expressions for the running of lepton mixing angles 
have been recently presented in Ref. [52]: When running M e g from the GUT scale down 
to low energies, the corrections to the leptonic mixing angle 0y are smaller than ~ \rrii + 
m j| 2 (l m «| 2 — I 772 ?! 2 )" 1 x where rrii and rnj are the eigenvalues of the ith and jth 

neutrino mass eigenstates of M e g- at the GUT scale. An appreciable running of leptonic 
mixing angles can thus only be expected in the IH or QD case. For NH neutrinos, however, 
the corrections are < 1° and, thus, negligible. Moreover, a tuning of phases always allows 
to switch off completely any RG effect on neutrino mixing angles - even in the case of 
inverse hierarchical and degenerate neutrinos [52]. A similar result has been obtained in 
the bottom- up approach in Ref. [53], where the starting point are the fixed low-energy 
observables. In addition, while the overall neutrino mass scale is affected by RG running, 
the neutrino mass ratios are hardly changed. Since our results should be very stable under 
RG running for the NH case (irrespective of the phases), we will focus in this paper on this 
type of hierarchy. 

4 Currently Allowed Realizations for the NH Case 

In this section, we focus on the constructed set of currently allowed seesaw realizations for a 
normal neutrino mass hierarchy (NH case). We obtain 173 084 different realizations for the 
NH case, which reduce to 8 030 cases if one does not count different phase combinations as 
different cases. This leads to 1 981 different texture sets, i.e., different combinations of Mg, 
M D , and M R . Naturally, we cannot show all of these possibilities in this paper. In Sec. 
we will therfore apply some selection criteria to reduce this dataset further and present the 
texture sets which seem to be most interesting to us. In this section, however, we discuss 
general features and some statistics of the constructed seesaw realizations for NH neutrinos. 

We concentrate on NH also because RG effects on neutrino mass ratios and mixing angles 
are expected to be small in this case (see Sec. 13.51) . This means that our generic assumptions, 
which may hold at some high energy scale (say at ~ 10 16 GeV), do not change significantly 
when running down to low energies where we match to experiment. In the NH case, one 
can easily diagonalize M e g for the allowed realizations in order to check that Eq. (1241) 
produces observables in agreement with current data. We have done this exercise: We 
have determined U u by diagonalizing M eff and computed £/pmns — U\U V employing the 
corresponding matrix Ug. From £7pmns> one can read off the mixing angles just as described 
in Ref. [25] o Note that we do not expect to reproduce exactly our input values in Eq. f)23p . 
since we do not require exact matching precision in Eq. ff24l) (see also Appendix IA. 21) . 

To describe the compatibility of a realization with current data, we use the performance 
12 See, e.g., also Ref. [54] for a discussion of neutrino mass matrix diagonalization. 
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Figure 4: Distribution of valid seesaw realizations for the NH case as a function of \ 2 as denned in 
Eq. (J2H). The values of #13 are all in agreement with current data. 
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Figure 5: Distributions of hierarchies in Mb (left) and Mr (right) leading to a NH neutrino masses for 
all valid seesaw realizations. 



2 _ /'sin 2 #12 - 0.3 V /sin 2 6 23 - 0.5 



indicator 

X = [ 0.3xcr 12 J + (, 0.5X6T 23 ) ' ^ 

(e.g., x 2 — 11-83 corresponds to a 3a CL exclusion for 2 d.o.f.). This corresponds to a 
Gaussian \ 2 approximation in sin 2 # 12 and sin 2 6*23 with the current best-fit values. For the 
relative la errors, we use a 12 — 9% (for sin 2 6 I 12 ) and <r 2 3 ~ 16% (for sin 2 # 23 ) [14]. Note that 
we only find sin 2 6 Vi 0.04 below the current bound, i.e., we do not have to impose an 
additional selection criterion. Fig. H] shows the distribution of the valid seesaw realizations 
as a function of x 2 defined in Eq. (1261) . Obviously, Eq. (124)) already ensures that the neutrino 
mixing angles of each realization are compatible with current bounds. It turns out that in 
all valid cases # 13 <^ 1° and only 6.5% of the realizations lead to 11.83 < \ 2 ^ 17 (which 
corresponds to a CL between 3 and 4a for 2 d.o.f.). Therefore, the selected realizations 
are all in perfect agreement with current data. Note that one might naively expect to find 
realizations with 613 around the best-fit value of 0° since this value has been used as an 
input for M^ p . However, in almost all cases one obtains # 23 to be around 50° despite of the 
best-fit input value of 45°. 

Fig. El shows the distribution of mass spectra or hierarchies proportional to (mf ,mf ,mf) 
and (mf,m^, m^) for and Mr respectively (each normalized to the corresponding 
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heaviest mass eigenvalue of Mp and Mr) for the NH case. These distributions are obtained 
by simply counting the number of realizations with a certain mass spectrum. Observe that 
M R has as a mass spectrum only (e 2 ,e, 1) and (e, 1, 1). Note that the spectra labeled as 
(e 2 ,e, 1) include also strongly hierarchical cases where the right-handed Majorana neutrino 
mass spectrum can actually be (e n , e, 1), with some suitable n > 3 (see Appendix lA.lj) . Fig. [5] 
shows that, for Mr, in all valid seesaw realizations, the mildly hierarchical mass spectrum 
(e, 1, 1) clearly dominates the (strongly) the hierarchical spectrum (e n , e, 1) with n > 2. As 
it will turn out later in Sec. this observation is also supported at the texture-level: more 
than 80% of the extracted textures lead to a mild mass hierarchy (e, 1, 1) for the right-handed 
neutrino masses. We will come back to this point in the next paragraph. The absence of 
a degenerate mass spectrum for Mr in the NH (and IH) case is a simple consequence and 
selection effect of Eq. f)27cl) in combination with the assumption M^ p ~ diag(e 2 , e, 1). There 
are many more possibilities for the mass spectra of M^ but they are dominated by the types 
(e, l,e) and (e, e, 1). In Fig. the pie piece "Other" also contains the hierarchy (e 2 ,e, 1), 
which implies that we have the same hierarchy in Mr (but not vice versa). In our method, 
no charged-lepton or quark-type hierarchy is produced. We find from Fig. [5] that there are 
many possibilities to obtain a normal neutrino mass hierarchy, but one cannot claim that 
this hierarchy appears typically in Mjj or Mr and then translates into M e g. 

The distributions of mass spectra in Fig. [5] may have immediate relevance for leptogenesis 
[55] when crudely extrapolating our results to the CP non-conserving case. In at least 80% 
of the cases that we found (c/., Sec. [5]), the right-handed neutrino mass spectrum is of the 
mildly hierarchical form (e, 1, 1). Thus, if the mass mf of the lightest right-handed neutrino 
is in the range mf < 10 12 GeV, the seesaw scale set by the mass of the heaviest right- 
handed neutrino would have to be significantly lower than the usual B — L breaking 
scale ~ 10 14 GeV. For the mildly hierarchical right-handed neutrino mass spectrum (e, 1, 1), 
successful leptogenesis might be achieved in two ways: (i) via resonant leptogenesis [56] 
(for recent models see, e.g., Ref. [57]) or (ii) by taking flavor effects into account [58] 
(for a connection with low-energy CP- violation see, e.g., Ref. [59]). In the resonant limit, 
mf could be as low as several TeV, thereby making this scenario testable at a collider. 
Strongly hierarchical right-handed neutrino masses, which is the standard case considered 
in the literature for leptogenesis, are in our analysis found to be by about a factor of 5 less 
abundant than the mild hierarchy. The possible strongly hierarchical right-handed neutrino 
mass spectra are all of the type (e n ,e, 1), where n > 2. Allowing n to be sufficiently large 
(say n — 8), the strongly hierarchical case can fit into a scheme with a seesaw scale of the 
order ~ 10 14 GeV and sufficient baryon asymmetry could again be generated through 
flavored leptogenesis. 

Fig. shows the distributions of mixing angles, where we concentrate on the number of 
maximal ("max") mixing angles = ir/A appearing in Up, Up/, Ur, and U#. If there 
is no maximal mixing angle, we call the scenario "All small". Let us first note that the 
distributions of the mixing angles in M v are very different from Ue. In Ug, we often find 
large mixings, which means that the large lepton mixing angles are not necessarily created 
in the neutrino sector, but can also come very often from the charged lepton sector. Note 
that the pie slice "All small" in Ue represents more or less CKM-like mixings in Ug. There 
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Figure 6: Distributions of mixings in Up, Ujj>, Ur, and Ue (in columns) leading to a normal neutrino 
mass hierarchy. The different pie labels refer to the number of maximal mixing angles, where "All small" 
corresponds to all mixing angles < e. 
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Figure 7: Fraction of special cases for normal mass hierarchy. The different pies show the fractions of 
the cases with symmetric Md, diagonal Mr, and Ue — 1, of all allowed realizations (not texture sets). This 
figure is based on the mixing matrices, where Ud — Ud 1 m the first case, Ur ~ 1 and Ue ~ 1 in the second 
and third case, respectively. For the similarity condition we allow e 2 -dcviations in the mixing angles. 
For instance, for an exact Ur = 1, one would have only 2% of all realizations. 

are also many possibilities with "trimaximal" mixinj^l (i.e., all three mixing angles 9fj are 
maximal for some given sector x) in Ue- This is different from the other mixing matrices, 
where trimaximal mixing hardly occurs, and either one maximal mixing angle or only small 
mixing angles are preferred. In particular, in Un> and Ur, only small mixings are typical. 

Fig. [7] shows the fraction of realizations that exhibit symmetric Mp, and/or diagonal Mr, 
and/or Ue ~ 1. In this figure, u Mp ~ symmetric" means that Md is symmetric up to 
possible corrections of the order e 2 . It is evident that there are only very few realizations 
with symmetric Mn or Ue — 1, and none with diagonal Mn, which is not surprising from 
what we have learned above. One may conclude from this result that there are plenty of 
possibilities to implement the seesaw mechanism without these constraints, which have, 
however, often been imposed in existing literature. 

Generally, note that, while one may argue that one can construct more possibilities in a 
more general and sophisticated framework, our realizations result from very generic and 
simple assumptions without adding another level of complexity. In this sense, our generic 
assumptions are much simpler than ad-hoc constraints, such as requiring that U& be diagonal. 

13 Not to be confused with tri-bimaximal mixing. 
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It is interesting to note that we find for Mr quite many realizations that are close to a 
diagonal form. 

5 A Selection of Textures for the NH Case 

In this section, we present a selection of texture sets for the NH case. These texture 
sets satisfy certain selection criteria listed in App. IB. II The full set of the texture sets is 
available in Ref. [60]. Table [T] shows these 72 texture sets for Mi, Md, and Mr, together 
with associated example realizations leading to these texture sets. For each texture set, we 
always choose the realization with the lowest x 2 , i.e., the realization which fits data best. 
For this realization, the table lists the mass spectra of Md and Mr, as well as the mixing 
angles 9fj in the different sectors. For each case in Table [U we have collected in Table [5] 
(in App. IB.2p all the corresponding phases (0 or it) in order to allowing for a complete 
reconstruction of the realizations and Yukawa coupling matrices. In addition, one can find 
there the PMNS mixing angles, as well as the number of realizations that become identified 
with each texture set through the texture reduction. In Table [H the parameter £ can take 
the values £ G {0, e 2 }. 
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Table 1: Complete set of selected seesaw textures/realizations for the NH case, where £ € {0,e 2 } (see text). 



In Table 12], we divide some of the interesting textures/realizations from Tabled] into certain 
classes, such as lopsided Mi (only #23 is maximal), anarchic0 Mp, (6^ = 9^ = #13 = #23 = 
j), lopsided Mp (at least one of the mixing angles 9f 2 and #23, x = D,D', is maximal - 
anarchic cases are excluded), hierarchical Mr (hierarchical mass spectrum but no maximal 
mixing ang le for M R , i.e., i< iag oc diag(e 2 ,e, 1)), semi-anarchic Mr (one of the angles 9^ 2 
or is maximal), "diamond" -type Mr = 7r/4 and the corresponding texture has a 
diamond shape [25]), and presence of a "dead angle £" (a mixing angle £ < e 2 that does 
not affect the corresponding matrix in the texture set) a Note that, as already mentioned 
in Sec. S] about 80% of all textures in Table [T] have a mildly hierarchical spectrum for the 
right-handed neutrino masses that is of the form M^ iag oc diag(e, 1,1), whereas a strongly 
hierarchical spectrum M^ iag oc diag(e n , e, 1) (n > 2) occurs in only roughly 20% of the cases. 
Although this classification may to a certain extent be incomplete, it could nevertheless serve 
to characterize significant features of the textures. Moreover, all realizations in Table [T] 
have in common that # 12 — 33°, #13 ~ 0°, and 6*23 — 51° (see Table Therefore, if future 
experiments measure 623 smaller than 45°, the presented realizations could be tested. 

Let us now illustrate how the textures in Table [1] could be generated in explicit models by 
considering two examples. For this purpose, assume an M-fold Zn product flavor symmetry 
group Gf = Illf^Zjy = Zn x Zn x • • ■ x Z^. We suppose that for each individual group 
Z$ there are two types of SM singlet flavon fields fa and f- that carry different Z$ charges 
but which are singlets under transformations of all the other groups Zff , with j ^ i. All 
flavons shall acquire universal VEVs: (/j) ~ (f-) ~ v, for i — 1,2, ... , M, where v — e ■ Mp 
and Mp is some fundamental Froggatt-Nielsen (see Sec. 12.21) messenger scale. Now, let us 
specialize to the case M = 7 and N = 4, and assume that each pair of flavons and // 
carries the Z± charges /j ~ 1 and f[ ~ 2, i.e., f t is singly and f[ doubly charged under 
zjfK Technically, this amounts to realizing fractional charges for the Z 2 -subgroups of Z± . 
We assign the leptons Gp quantum numbers as shown in Table [3] for two example models. 
In both models, we have assigned in Table [3] each lepton a row vector (qi, qi, . . . , qi), where 
qi denotes the Z 4 charge of the lepton under the group Z± (i = 1,2, ... ,7). Models 1 and 
2 in Table [3] respectively produce the texture sets #17 and 18 in Table [1] via the Froggatt- 
Nielsen mechanism. Since we have already found in Tables [1] and [5] valid realizations for 
these textures, we know, without any further calculation, that the order unity Yukawa 
couplings in model 1 and model 2 can be chosen such that they reproduce the lepton mass 

14 Cases of anarchic Mr cannot appear in Tabled] due to the selection criteria in App. IB. II 
15 For example, if a matrix element reads A e 2 +B £e+C(e 4 ), where A and B are order one coefficients, then 
the choice of £ £ {0, e 2 } has no impact on the texture, i.e., A e 2 + C(e 4 ) — > e 2 and A e 2 + B e 3 + C(e 4 ) — > e 2 - 
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Class Texture # 



U t * Vckm 46, 47 

U t ~l 33, 34, 71, 72 

0f ^ f 20, 27, 33-35, 44, 46, 47, 49-53, 57, 58, 62-64, 

67, 68, 71, 72 

Lopsided M t 1, 4-8, 10, 11, 13, 43, 48, 56, 60, 61, 69, 70 



BimaximalM^ 
Trimaximal Mi 


9, 17, 21-23, 26, 32, 38-40, 55, 59, 65, 
3, 12, 14, 15, 16, 18, 19, 25, 28 


U D - u D > 

Anarchic Mo 
Lopsided M D 


12, 23 
24, 54 

8, 20, 21, 23, 27, 32, 33, 35, 41, 42, 45, 49, 
50, 52, 53, 62-64, 67, 68 


Hierarchical Mr 
Semi-anarchic Mr 
Diamond Mr 


10, 33 

48, 58, 65, 66 

1, 16, 18, 24, 32, 46, 47, 49, 50, 54-57, 60-63 

2, 3, 5, 9, 15, 17, 19, 29-31, 35, 43, 44, 67-72 


ilC g oc M* ag 
i< iag oc diag(e 2 , e, 1) 


43, 44, 48, 58, 65, 66, 71, 72 

1, 43, 44, 48, 49, 50, 58, 61, 62, 63, 65, 66, 71, 72 


Dead angle £ 


1, 6, 7, 9, 10, 13, 14, 24, 28, 30-33, 36-38, 41, 42, 
46-48, 52, 53, 57, 61-64, 68, 70-72 



Table 2: Classification of textures and realizations which may be of special interest. Here, "~" means 
up to phases. The case M^ lag oc M^. lag allows, in combination with Eq. (|27c[) . only for the mass spectra 
(e 2 ,e,l). Note that the relative fraction of hierarchical M^, lag in Tabic Q] is much higher than prior to 
applying the selection criteria, as it is evident from Fig. [5] 



and mixing parameters in perfect agreement with data. Indeed, the explicit realizations in 
Tables [T] and [5] allow, if one wishes, for a complete reconstruction of such a valid set of order 
one Yukawa couplings. 

The two models above represent only very specific examples of how one could directly 
apply Table [Tj to identify the possible flavor symmetries and their breaking in a model. 
However, there are certainly many more possibilities. For example, we assumed here, for 
simplicity, copies of a single discrete Abelian flavor symmetry group, but new possibilities 
arise when considering product groups with different Zn subgroups or also non-Abelian 
flavor symmetries. Moreover, note that, since U# has been set equal to the unit matrix, 
Table[TJlists only a few percent of the actual total number of texture sets in extended QLC. A 
systematized scan for models generating the textures in extended QLC should also include 
these extra cases. In addition, it would be interesting to explore the compatibility with 
further constraints from GUT-relations or cancellation of anomalies (see, e.g., Ref. [61]). 
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Field 


Model 1 


Model 2 




(0,0,0,1,0,1,1) 


(2,0,0,2,0,0,1) 




(2,0,0,1,0,1,1) 


(2,0,0,2,0,0,1) 




(0,0,0,1,0,1,1) 


(0,2,0,0,2,1,0) 


f) 


(0,2,0,0,1,0,1) 


(2,0,2,2,2,1,0) 




(0,0,0,0,1,1,0) 


(2,2,0,2,2,1,0) 


h 


(0,0,2,1,0,0,1) 


(0,2,2,2,2,0,1) 


e\ 


(2,2,2,1,1,1,1) 


(0,0,0,0,0,1,1) 




(2,2,2,0,0,1,0) 


(2,2,2,0,0,3,3) 


4 


(0,2,2,0,0,0,0) 





Table 3: Assignment of Gf charges to the leptons. Models 1 and 2 lead respectively to the textures #17 
and 18 in Table Q] (see text). 

6 IH and QD Case 

Even though the main focus of our study is the case of NH neutrino masses, we have also 
calculated the valid realizations for IH and QD neutrinos using Eq. f|T2|) for the respective 
mass hierarchies. These cases are qualitatively different from the NH case because RG 
effects may be relevant here. In this section, we present a qualitative discussion of the IH 
and QD cases. 

6.1 IH Spectrum 

In the IH case, we find 797928 different realizations, which corresponds to 10 196 different 
qualitative cases ignoring phases. This is about a factor of five more realizations than in 
the NH case. After texture reduction, one obtains 4 268 different texture sets, i.e., different 
combinations of Mi, Mp, and Mr, which is about twice the number as in the NH case. 
This indicates a higher redundancy at the level of the realizations as compared to the NH 
case, which we will comment on more explicitely in the QD case in Sec. 16.21 Since some 
of the selection criteria for NH neutrino masses are based on the \ 2 selector, which is not 
defined in the IH case, the texture reduction of the texture sets for the IH case has to be 
carried out in a way that differs from the procedure described in Appendix IB. 11 However, 
by arguments similar to those in Appendix IB. 11 we find a reduction of the number of cases 
that is comparable to that in the NH case. 

The mass spectra of Mr> and Mr for the IH case are shown in Fig. [SJ In contrast to the 
case of NH neutrino masses (c/., Fig. [5]), Mr, shows a larger variation of possible mass 
hierarchies. It turns out that most realizations for Mr, have a mass hierarchy of the type 
(e, e, 1). In Mr, we obtain the same mass spectra as for the NH case, but with a nearly 
opposite weighing. For Mp, one may have an inverse mass hierarchy, whereas M^ iag has by 
construction a "normal" ordering (c/., Appendix [ATT]) . 

Fig. [9] shows the distribution of maximal mixing angles for the IH case. Note that one 
observes a similar distribution for NH neutrinos (cf. Fig. [6]), namely that in most cases large 
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Hierarchies in M D 



Hierarchies in M R 




Figure 8: Distributions of mass hierarchies for Md (left) and Mr (right) leading to IH neutrino masses 
for all seesaw realizations. 



Figure 9: Distributions of mixing angles in Ud, U' d , Ur, and Ui (in columns) that lead to IH neutrino 
masses for all seesaw realizations. The different pie labels refer to the number of maximal mixing angles, 
where "All small" corresponds to all mixing angles < e. 

leptonic mixing angles come from Ug. However, for IH neutrino masses, there is a slight 
tendency to larger mixing angles in Ud, Ud 1 , and Ur, and smaller mixing angles in Ui than 
in the NH case. The distribution of the special cases in Fig. [7] resembles that for the normal 
hierarchy. For example, Mr is approximately diagonal in about 15% of all cases. 

6.2 QD Spectrum 

The case of QD neutrino masses is more tricky from the computational point of view because 
of redundancy. As a simple example, consider Eq. (120bl) with M eff ag = 1. In the real case, 
one can use any C/pmns and Ue in this equation in order to obtain M^ p = 1, i.e., the different 
cases for Ue are not really qualitatively different and a statistical analysis, such as above, 
does not make much sense. In addition, the matching to experiment (or the experimental 
data as a selection criterion) is, in this way, rather meaningless. From Eq. ( liHl) . one can read 
off that we require oc l + 0(e 3 ). This implies that any realization producing ~ c- 1 
will be accepted by the algorithm, which are quite many. An interesting observation in the 
degenerate case, is that the realizations divide into two qualitatively different cases: the 
trivial case, i.e., Md = Mr = 1, and a large number of non-trivial cases. There is no such 
trivial case for IH neutrinos. 



Mixing angles for inverted hierarchy 
U D U d U r 
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7 Summary and Conclusions 



Quark-lepton unification generally leads to predictions in the patterns of the Yukawa cou- 
pling matrices of the fermions, which are often called "textures". Because of the high 
dimension of the parameter space, a direct systematic study or scan of all valid textures at 
the level of the Yukawa couplings is difficult to perform. Therefore, we suggest a "bottom 
up" search strategy: We construct the parameter space of all valid textures from very generic 
(low-level) assumptions in the relevant leptonic mixing matrices and mass spectra. In this 
approach, we use the context of extended quark-lepton complementarity, which means that 
all mixing angles are either maximal or described by powers of e ~ 9 c- In addition, all mass 
ratios (of the neutrinos and charged leptons, as well as the mass hierarchies in Mn and Mr) 
are also parameterized by powers e. This expansion parameter may be the remnant of a fla- 
vor symmetry describing the masses and mixing angles in both the quark and lepton sectors 
within a quark-lepton unified context. Note that in this setting, the solar neutrino mixing 
angle can only emerge as a combination of maximal mixing and 9c, i.e., our assumptions 
are somewhat more general than the usual QLC relations [see, e.g., Eq. (jSJ)]. 

Our procedure, in short terms, is as follows (cf., Fig. [2]): First, we systematically construct 
all combinatorial possibilities for all mixing matrices and mass spectra for the type-I seesaw 
mechanism up to the order e 2 . In each case, we compute the effective neutrino mass matrix 
from that. In the second step, we determine all possibilities for M^ p containing the 
experimental information from £/pmns ; the neutrino mass hierarchy, and Ug. In the third 
step, the two matrices are matched, i.e., the realizations compatible with current data are 
selected. Finally, we identify the leading order entries in the matrices, which lead to the sets 
of textures. We scan about 20 trillion possibilities focusing on the real, CP conserving case, 
and we mainly discuss a normal neutrino mass hierarchy. As an experimentally motivated 
input value for sin 2 2# 13 , we take the current best-fit value sin 2 2# 13 = as an example. 

Compared to individual models in the literature, our assumptions are too generic to be 
able to predict the outcome a priori, which may reduce the bias. This key feature allows 
for the interpretation of the valid result statistics with respect to observables, constructed 
hierarchies, etc.. For example, we find many different cases with a very mild hierarchy in 
Mr (with an abundance of about 80% in the textures), a result which does not involve any 
bias from the point of view of the assumptions. In addition, we find only very few cases 
with a symmetric Dirac neutrino mass matrix M^ or small mixing (Ug ~ 1) in the charged 
lepton sector. However, we find a roughly diagonal right-handed neutrino mass matrix Mr 
in relatively many cases, i.e., in about 24% of all realizations for the normal neutrino mass 
hierarchy. It would be interesting to investigate signals of lepton flavor violation for our list 
of matrices, and connect our results to statistical studies along the lines of Ref. [62] or to 
recent attempts of scanning the SM parameters including quark and lepton masses [63]. 

We have presented a complete selection of 72 examples for Yukawa coupling textures sat- 
isfying specific selection criteria. A more complete list of 1981 texture sets (obtained by 
relaxing the selection criteria) can be found in Ref. [60]. We have shown in two examples 
that the textures are very useful for a direct search for flavor symmetries predicting current 
data. Therefore, our list(s) of textures can be understood as an intermediate result from 
the point of view of model building. They might be used to systematically test and find 
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flavor symmetries and allow for a more systematic approach to constructing models. 

We conclude that systematic, machine-supported searches of parameter spaces related to 
model building may, in fact, open up new possibilities. While conventional approaches focus 
on individual models in greater depth, our approach produces all valid lepton mass matrices 
with the only bias of the (often very generic) input assumptions. Therefore, they allow for 
the identification of new possibilities, as well as for more general studies of the discussed 
parameter space. Hence, our study should not only be interpreted with respect to the input 
assumptions, but also with respect to the procedure itself. 
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A Details of the Procedure, Heuristics, and Complexity 

In this appendix, we give details of the algorithm, heuristics (also physically relevant ones, 
everything which makes the algorithm faster but does not lead to a loss of generality), and 
complexity. 

A.l Parameter Space for the Mass Hierarchies 

In the procedure described in Sec. 13.3} we have considered all possible diagonal mass matri- 
ces M^ iag = m D diag(e a , e b , e c ) and M^ iag = m D diag(e a ', e b ' , e c '), with suitable non-negative 
integers a, b, c, a', b', and c', that are compatible with the observed neutrino mass squared 
differences and leptonic mixing angles. As we will now see, the relevant range of these inte- 
gers can actually be restricted considerably if we take the parameterization of the neutrino 
mass spectrum in Eq. (fT2"1) into account. 

By factoring out common powers in e, we see that M D ia,g and M^ iag can be written as 

M* ag = m D diag(e"\ e n , 1), Af£ iag = m D diag(e m , 1, e n ), M* ag = m D diag(l, e m , e n ), 

(27a) 

and 

M^ iag = Mb-l diag(e p , e\ 1) with < q < p, (27b) 

where m,n,p, and q, are non- negative integers, and mo ~ 10 2 GeV. In Eq. (I27bl) . we 
have made use of the possibility to bring M^ iag to a strictly hierarchical form, which is 
expressed by the condition < q < p. One then has, however, no longer the freedom to 
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choose the order of the mass eigenvalues of M D mg , and that is why Eq. (|27aj) includes all 
permutations of M D iag without any specific ordering of m and n, i.e., we allow in Eq. (I27al) 
for both the cases n < m and n > m. In the following, we set m v = m 2 D / Mb-l- Although 
mo ~ 10 2 GeV and Mb-l ~ 10 14 GeV would be preferred values, we can, of course, always 
rescale mo — > a ■ nip and Mb-l a 2 • Mr by some factor a, thereby leaving the absolute 
neutrino mass scale m u ~ 1CT 2 eV unchanged. From detM^ ag = (det M^ iag ) 2 det (M^ iag ) _1 
we find the relation 

2(m + n) — p — q = k, (27c) 

where k is a non-negative integer that depends on the type of neutrino mass spectrum in 
Eq. ffl2|) : We have k = 3 for a normal, k = 1 for an inverted, and k = for a degenerate 
neutrino mass spectrum J^l 

We are now interested in determining the allowed ranges for the integers m, n, p, and q, 
in Eqs. fl27|) . Later, we will apply Eq. (124"I) up to order e 2 , i.e., we will require a nu- 
merical matching precision 0(e 3 ) between and M^ p . From that matching preci- 
sion it follows that only powers up to e 2 are relevant in the individual factors in the 
product Eq. fl20al) . because higher order terms will be absorbed by this matching uncer- 
tainty. The only non-trivial aspect in this argument is the factor M^ 1 . Let us write 
(M^ iag ) _1 = e- p M~ l _ L diag(l, e p ~ q , e p ), and redefine m D -> e p ' 2 m D and M R -> e v M R , which 
leaves m u unaffected. Consider now the entries e m , e n , e p ~ q , and e p , in M D iag and (M^ iag ) _1 . 
If m, n, p — q, or g, are larger than 2, then the corresponding contribution to M e R will be 
absorbed in the matching precision. It follows that it is sufficient to restrict the maximum 
values of m,m,p, and q to 2, i.e., we need to consider for the powers of e only 

< m,n,p,q < 2, (28) 

where we have used that < q < p. The matrices M D 13lg and M]^ iag with higher powers 
of e fall then into one of the classes already covered by a combination of powers satisfying 
Eq. (J22|. 



A. 2 Matching of Matrices and the Absolute Neutrino Mass Scale 

In this section, we describe the numerical implementation of Eq. ($M§: M^\ t =o.2 — M^ p \ e=Q ,2- 
For our procedure, it makes sense to require a matching precision of 0(e 3 ) since e 2 is the 
maximum power that is used for the mass hierarchies and mixing angles (which was chosen 
because higher orders would be absorbed by the current measurement precision). This im- 
plies that a higher order matching precision will be too precise/restrictive because e 3 powers 
are not directly produced, and a lower order matching precision would not be able to dis- 
tinguish cases differing by e 2 -terms. Numerically, the 0(e 3 ) term can have any order one 
coefficient. We require a precision of P = e 3 for Eq. (1241) . which turns out to be reasonably 

16 From Eq. (|27c[) and p > q, it also follows that m,n < | + p. 
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restrictive by using actual tests, i.e., we require 17 ! 

(M^ii - P < c ■ (AfS)^ < (M e 7) 4 , + P, (29a) 

with P = e 3 and c is a constant which can be absorbed into (or can come from) the absolute 
mass scale, mixing angles, and the highest hierarchy power in Mr. 

In order to evaluate Eq. (129a[) . we have to determine c in order to check whether falls 
into the proper interval. Let us re-write Eq. (I29al) as 

(M ex P) low < c . (M th) . . < (M ex P) up (2%) 

This corresponds to 6x2 independent inequalities. Dividing now these inequalities by 



W|< c <wi for(M , h , >0 



C /If exp ) np ( A/f cxp "1 low 

1 eff iM < r < cff Ai for ( M th ) < (2Qc) 

Watch for the special case (M e g)ij = 0. The intersection of these intervals is then obtained 

as 

max I nTrth^ ) < c < mm I ( J I (29d) 

S v ' 

with (M^ p )^ and (M^ p )^- chosen according to the case selection in Eq. f l29cf) for i and j 
individually. Only, if c L < c u , then we have an allowed range for c. If, however, c L > c u , 
the realization is refused. Note that the constant c is not used anymore further on since it 
may have different origins. However, the appropriate absolute mass scale will be implicitly 
determined in order to satisfy Eq. ( T24l) . 




A. 3 Complexity and Counting of Cases 



In order to avoid in our procedure the generation of equivalent realizations and double- 
counting of cases, several heuristics can be used. As far as the phases in Eq. f)20ap are 
concerned, we only have the matrices Dp, K, D for the real case, which leads to 2^ 3+2+3 ^ = 
2 s = 256 phases. For 6f 3 = 0, the phase 8 X is unphysical, and we only test one case. This 
leads to 4 6 x (3 x 2 + l) 3 different angle combinations in M% For M^ p in Eq. f!20bl) . we 

17 The dataset of valid models is very sensitive to this P. If it is much smaller, we do not find any valid 
realizations because e 3 -terms are not directly produced by our procedure. If it is somewhat larger, the 
found realizations fail consistency with data. For example, for a normal neutrino mass hierarchy, one can 
diagonalize and read off the mixing angles as in Ref. [25] . These mixing angles can be compared with 
the actual ones which have been used as input values for M^ p . We find that our matching precision is 
sufficient such that only realizations compatible with current data at least at the ha confidence level arc 
allowed for the NH case (most of them actually fit much better); c/., Fig. 0] A weaker constraint on P 
allows more realizations which provide, however, a too bad fit. 
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do not generate Dg in Eq. (120bl) . because the phases in the corresponding Dp will match 
a valid case of M^ p since both matrices appear on the outside (a valid phase in Dp then 
actually corresponds to two cases for Dp and Dg). In the real case, we therefore have only 
2 2 = 4 phases from Kg, and 4 2 x (3 x 2 + 1) = 112 cases for the angles in Ug. This has to be 
multiplied with the corresponding allowed hierarchies in Mr, and Mr. It is easy to derive 
that there are 12 possible cases for a normal neutrino mass hierarchy, 9 possible cases for an 
inverted hierarchy, and 13 cases for the degenerate case (c/., Appendix IA.1I) . In addition, 
we test six cases for the true values in {Jpmns- Since in the degenerate case any Ug and any 
fpMNS leads to a diagonal M^ p , it is sufficient to test one case for Ug and Kg. This means 
that the degenerate case hardly contributes to the complexity. Further simplifications are 
possible for M^ iag = 1, etc.. In this case, one has to watch that the relative counting of 
different hierarchy /mixing angle cases is affected. In total, we test 

N = 256 x 4 6 x (3 x 2 + l) 3 x (112 x 4 x (12 + 9) + 13) x 6 ~ 20 • 10 12 (30) 

different combinatorial possibilities, which require about 2 months of running time on a 
modern computer using a C-based software. 

The number of allowed combinations may be used as a statistical measure in some cases. 
The interpretation is then related to the number of different possibilities. 

B Details of the Results for the NH Case 

This section contains supplementary material for Sec. EJ such as the selection criteria used 
and the phases necessary for a complete reconstruction of the charged lepton and neutrino 
Yukawa couplings. 

B.l Texture Set Selection Criteria 

Let us now describe the selection criteria leading to the texture sets presented in Sec. |5j In 
producing Table [U we have used the following criteria: 

1. The seesaw realizations should resist an increased experimental pressure provided that 
the current best-fit values are unchanged, i.e., sin 2 2^13 will not be found. We impose 
an extrapolated experimental limit by using a\2 — 4.6% [44] and 023 — 10% [43] in 
Eq. (1261) . which corresponds to about a decade from now. Since # 13 is in any case 
smaller than 1° and would anyway resist an increased experimental pressure, we do 
not have to consider it for a further selection. 

2. Textures that differ by entries of the order e 2 should be included with only one example. 
We choose the realization with the best \ 2 ■ 

3. The selection of realizations should be stable under e 2 - variations of the mixing angles, 
i.e., if we find mixing angles e 2 , the realization has to be valid for e 2 — * as well. We 
include the newly obtained texture sets in Tabled] by introducing the angle £ G {0, e 2 }. 
Note that these cases are initially generated as well, but they might be filtered out 
later by Eq. 
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Selection criterion 


^Textures 








None 


1 981 


20 


621 


35 


1.-2. 


1 048 


20 


475 


31 


1.-3. 


447 


20 


270 


22 


1.-5. 


72 


17 


65 


21 



Table 4: Total number of texture sets, and the numbers of distinct textures for Me, Md and Mr. Each 
row corresponds to the application of the selection criteria specified in the first column (see text). 

4. We omit texture sets with anarchical Mr (matrices just filled with entries "1"), since 
such structure-less textures do, in general, not yield much useful information. 

5. We also show only one example for each texture set (Mi,M 2 ,M 3 ) ~ (Mi, Md, Mr), 
where Mi (i = 1, 2, 3) are the corresponding textures obtained after texture reduction, 
which have one pair {Mj,Mk} appearing together with more than one possible Mi 
(i ^ j ^ fc). In this case, we keep the texture set associated with a representation 
that has the lowest x 2 ■ 

By using the above selection criteria, we obtain the list of 72 texture sets shown in Tables [1] 
and [5j The effects of the different selection criteria on the number of texture sets is shown 
in Table HI In this table, we also show the numbers of distinct textures for Me, Mo and 

Mr. 

B.2 Supplementary Information for the Texture Sets 

Table [5J shows extra details of the realizations listed in Table [TJ It contains the complete 
set of phases, the PMNS mixing angles, x 2 (i n the 10 years limit, cf., App. IB. II) . and the 
number of realizations leading to each texture set (for £ = e 2 in ambiguous cases). Note 
that we chose iff = ip® = ip® = af = a® = since these phases appear in M^ only in 
combination with other phases (see Eq. (j20p ). and can be absorbed into the other phases. 
Tables [T] and [5J provide together the complete information sufficient to fully reconstruct the 
Yukawa coupling matrices of the 72 realizations. 
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Table 5: Supplementary information for Table Q] (fi'2 3 = Oi\2 = 0) ( see text). 
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